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1 Introduction 



We investigate the critical behaviour of hermitean matrix models in the double 
scaling limit N — > oo and g — ► where N x N is the size of the matrices, 
g is a coupling parameter, and a finite number / of such random matrices is 
coupled to a chain. This scaling behaviour describes two-dimensional quantum 
gravity coupled to the matter fields of rational conformal field theories For 
/ = 1 and / = 2 the existing analysis of double scaling behaviour is complete 
P], For / = 3 only a small number of examples are known leading to just 
three universality classes. We are going to present an algorithm which allows the 
systematic construction of all double scaling limits for / > 3. For an infinite 
subclass of / = 3 models this algorithm employs only linear algebra and all 
critical coefficients are rational. For all other cases of / = 3 and in particular 
for / > 4 we need to solve zeros of polynomials and the critical coefficients turn 
into algebraic numbers. We emphasize that matrix models are deeply related 
with Toda hierarchies 0] but that this relationship has not been fruitful for the 
elucidation of the critical behaviour. 

In order to fix the notations we briefly describe hermitean matrix models. The 
action is 

r / 

Ct=l 

-J2c a M ia) M {a+1) \ (1.1) 

a=l ' 

where each potential V a is a polynomial of degree l a (only l a > 3 is of interest). 

W) = X>K (L2) 
fc=i K 

Throughout the paper we assume that in relations between l x and l 3 w. 1. o. 
g. h > I3, if not stated otherwise. Stability will be completely neglected in this 
work. The partition function is 

r f 

Z= dM^e~ s (1.3) 

dMM = J] d(ReMlf)d(ImM { k ?) (1.4) 

i < j 
k < I 
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The method of orthogonal polynomials || |9] makes use of biorthonomal systems 
of polynomials 

{n m a , n m // } 

A,/i6E (1.5) 



satisfying 



x 



f[ dA^n m (A«)n n (A^) 

exp { - £ V a (A< a >) + £ c Q A^A( a+1 ) 

q=1 a=l 



= ^mn (1-6) 

Differentiation and multiplication matrices are introduced by 

n 

K = J2( A f)n m K n (1.7) 

n 

AIT m (A) = E( 5 0mnn n (A) 

n 

fjJI m (iJ,) = ^2(Bf) nm fl n (ii) (1.8) 

n 

With the help of auxiliary matrices 

B 2 , B 3 , ...Bf_2 

we can derive Dyson-Schwinger equations 

A 1 + c 1 B 2 = Vl(B 1 ) 
c a-iB a -\ + c a B a+ i = V^(B a ) 
2 < a < f - 1 
Af + c f - 1 B,- 1 = V}{B f ) (1.9) 

These matrices have support at 

/ 

{A\) mn = except for — ]^[ (/q, — 1) < n — m < — 1 

a=l 

/ 

(Af) mn = except for 1 < n - m < JJ (l a - 1) (1.10) 



and 



(B a ) mn = except for - J[ (lp - 1) < n - m < J[ (l p - 1) 

0>a f3<a (1.11) 



From (O), (Ofi follows 



[B 1 ,A x ]=l (1.12) 

and the Dyson-Schwinger equations imply 

[B U A X ] = ci [B 2 ,B X ] =c 2 [B 3 ,B 2 ) = ... 

= c f -i[B f ,B f _ 1 ] = [A f ,B f ] (1.13) 



One can easily scale all the {c a } to one in (|1.9|), ( |1.13| ) 



For convenience we present a summary of our results at the end of this work. 
However, we can anticipate that the universality classes found are those and only 
those [p, q] of the two-matrix model: p and q are either coprime or coprime after 
a division by a factor r different from p and q . Then p and q are the orders of 
a pair of differential operators of a generalized Korteweg-de Vries hierarchy. 



2 Solving the Dyson-Schwinger equations 

The solutions of the Dyson-Schwinger equations are obtained in two steps pertur- 
batively. One observes that any singularity obtained in the double scaling limit is 
already determined by the solution of the Dyson-Schwinger equations at leading 
perturbative order. At this order we assume 



(B a ) n , n+m — p$ (2.1) 
l.o. 

so that for each B a we obtain a generating function 

( a )( -y\ — \ " ni )?™ 1 



r v ,< z 



E/i^ m (2.2) 



with z a complex variable (see ( |2.17| )). The are submitted to the constraints 



( p.. 11] ) so that r^ a '(z) is rational with poles only at z = and z = oo. These 
generating functions are inserted into the "internal" Schwinger- Dyson equations 
(/>3) ' 

B a -i + B a+1 = V^(B a ), 2 < a < f - 1 (2.3) 
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Each r^ a \z) is required to exhibit a zero of order X a at z = 1. The ansatz that 
leads to a set of zeros {A Q } is called a "maximal critical point" if it fixes all 
unknowns in the Schwinger- Dyson equations at leading order (no parameters are 
left over). We determine only maximal critical points. The {pff} an d the critical 
coupling constants (|1.2p 

{g^}, 2<a</-l (2.4) 

can then be calculated. In the subsequent sections 3 (for / = 3) and 4 (for / > 4) 
this program is performed. The {pff} are either rational or algebraic irrational. 
In the case / = 3 there is a whole subclass of solutions where this program can 
be executed analytically resulting in only rational solutions. 

In the second step we consider the "external" Schwinger-Dyson equations 

A 1 + B 2 = V{{B X ) 
A f + B f ^ = V}{B f ) (2.5) 

They are (at this perturbative order) only used to fix the remaining critical cou- 
pling constants 

{^ a) }, a = l or a = f (2.6) 

This is done as follows. Denote the restriction of to its holomorphic part 

at zero (infinity) by 

r^(z)>, (r^(z)<) 
respectively. Then using (|1.10| ) we can reformulate (|2.5|) as 

r®(z)> = V{^ l \z))> 
r^- x \z)< = V' f (r {f \z))< (2.7) 



Now from ( 1.11]) 



r W(z) = p^z + r (l \z)< 

r V\z) = p { l\z- 1 + rV\z)> (2.8) 
(pS x) ^ 0, p U \ ± 0). 



We see that we can introduce new variables 



f\ = r^ l \z) 



f-J 1 = r u \z) (2.9) 
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by holomorphic maps in a neighborhood of z = oo (z = 0), respectively, conse- 
quently 

r^(z)i = fl + Oi^r) (fi->oo) 

r</>(*)* = fj k + 0(f f ) (f/-0) (2.10) 
and by a series expansion we get 

n 

fc=— oo 
+00 

^ n = E a 2?(*v) fc » ^<o (2.11) 



k=n 



The {a^} are rational functions of the {p^} and can be calculated recursively. 
Then the desired critical coupling constants result from 

o (1) - YVV 1} 

m=0 

<7& = E P^^S (2.12) 

m=-(/ / -l) 

Having performed the leading order solution, the perturbative expansion pro- 
ceeds as in the case / = 2. We keep the critical coupling constants fixed, multiply 
the whole action with 

N 
9 

and tune 

N^oo, g->g c (2.13) 

as follows. 

The matrix labels n, m become continuous in this "double scaling" limit 

Ti 

jj = t 0<e<l (2.14) 
We replace the label N by the "string coupling constant" a : 

l = a 2 ^, 7 <0 (2.15) 



so that a — > for N — > 00. Moreover 
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t= 9 -{l-c?x) (2.16) 



The variable z ( |2.2| ) is dual in the Fourier series sense to the discrete matrix label 
m (see (|2.1|)). Now we set 

z = e iLp (2.17) 
<p = a~~<p (2.18) 



so that due to ([H|), ftTHD 



is the quantum mechanical momentum operator corresponding to x. The pertur- 
bative expansion is in powers of 

and 7 is defined at the end from the Korteweg-de Vries hierarchy operators by 

7 = ~ 2 (2.20) 
P + g - 1 

In this limit the matrices {-B Q }^ =1 become differential operators {Ra]^ a=l 

B a ^a~^R a ae{l,...,f} (2.21) 
where the order of R a is \ a . While the /-matrix model belongs to a class denoted 

(Zi, l 2 , ■ ■ ■ , If) 

we use square brackets 

[Ai, A 2 , . . • , A/] 

to denote the critical points. 

3 The three-matrix model 

For the three-matrix model we proceed as follows. We set 

(-Bi) n>n+fe = r^\n) 
(£ 2 )n,n+fc = rf\n) 

(B 3 ) n;n+k = rf (n) (3.1) 
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and make the ansatz 



r£\n--) = p^+a-^ix) 
rf\n-\) = pf+a^ufix) 

rf(n-i) = pf + a-^uf{x) (3.2) 

By inserting this ansatz in the Dyson-Schwinger equations we get a perturbed 
system in powers of a -7 . 

To 0-th order one has equations between the p^\ p k 2 \ P^ an d the coupling 
constants g^\g^\g^'. First we will concentrate on this system of equations. 

In order to solve them we make a more precise ansatz for the pj^ , p^ and pj^ . 
We introduce the generating functions 

r«(s) = E p^z m 



m=-(h-l)(l 3 -l) 
Pi-1) 

r (2) (-) = E 

m=-(Z 3 -l) 
(l 3 -l)(ll-l) 

(3)/^ _ \^ „(3)_m 



r 



^) = E ^ (3-3) 



-l 



and assume they are of the form 



r m iz) 



(z- l) Al P fe ^i )fe _i) + i_ Al (^) 

(z-l) x *P {ll+l3 - 2 y X2 (z) 



z 



(Js-l) 



r (3) (z) = (^-l) A3 ^i-l)fe-D+l-A3^) (34) 

where the P n are polynomials of degree n. 

Inserting such an ansatz in the Dyson-Schwinger equations of order zero one 
obtains the values for the critical coupling constants gk)9k>9k anc ^ can ^ x ^ ne 
coefficients of the polynomials P n . 

In the case of the three-matrix model, in contrast to the two-matrix model 
(compare 0), more than one maximal critical point can be found for each model. 
For example the (4,3,3) three-matrix model has three different maximal critical 
points. They will be discussed in detail in Appendix B. 

Among the critical points exist in general two different types which will be 
called type I and type II. Critical points of type I are given if Ai = A3. Then we 
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have R\ = —R 3 and the commutators [R 2 , R\] and [R3, R2] are obviously identical 
and therefore compatible with (|1.13|) . For the critical points of type II we have 
Ai 7^ A 3 . If Ai < A 3 then to leading order Ri and R 2 are identical or R 1 is a power 
of R 2 up to a multiplicative constant and their commutator is zero. But in higher 
orders the commutator of R\ and R2 is the same as the one between R 2 and -R3 
and again [R 2 , Ri] and [R3, R 2 ] are compatible. Of course, in this case the role of 
Ri and R3 can be exchanged. 

The appearance of the two types can be understood from examining the "in- 
ternal" Dyson-Schwinger equation (|2.3| ). We may choose the ansatz such that to 
leading order Bi and B 3 compensate each other (type I) and V 2 '(B 2 ) only con- 
tributes in higher orders. Or to leading order B\ and V 2 (B 2 ) compensate and B 3 
only is important for higher order contributions (type II). 

In particular with the ansatz 

\ (2 - l)^ 3 " 2 ) , . 

^ (2) W = K ~^h) ^ 

which yields maximal order for the differential operator R 2 , and for which the 
Dyson-Schwinger equations can be solved by linear algebra only (see Appendix 
A), we obtain several interesting cases. If (Zi — 1) and (I3 — 1) have no common 
divisor then we found a critical point of type [l 2 + 1, l\ + l s — 2, l 2 + 1). In the 
other cases, in which n{l\ — 1) = m(l 3 — 1), n,m G IN\{0} holds, we got critical 
points of higher order than [l 2 + 1, l x + 1 3 — 2, l 2 + 1] except in the cases where 
(li — 1) = m(l 3 — 1), m G IN\{0, 1} and l 2 > m. In these cases the construction of 
the maximal critical point fails which is shown in detail in Appendix A. 

In the following table we have listed some three-matrix models and their max- 
imal critical points. Rational solutions are marked with the abbreviation "rat". 
Moreover we give as last entries in this table the universality class [p, q] of these 
critical points which were defined first for two-matrix models. 
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L, l>2, h) 


[Ai, A2, A3J 
nist type 


r \ \ \ 1 

|Ai, A2, A3J 

second type 


fa, 


/( o q 


/I /I T.I t" 

4,0,4 rax. 




c; /i 
0,4 






A A K 

4,4,5 


5,4 






Q Q P. 
0,0,0 


7 Q 
' !<3 


A A 1 


5,5,5 rat. 




ft 17 
0,0 




K Q ft 
0,0,0 




7 Q 
' ,0 






A A P. 

4,4,0 


P. A 

0,4 






Q Q 7 

3,3,7 


7 Q 

7,3 


/I i; o 


6,5,6 rat. 




ft 17 
0,0 






4,4/ 


7 /I 






Q Q 8 
0,0,0 


8 Q 
0,0 


A ^ 


1,0,1, 1 at . 




7 ^ 






4,4,8 


9,4 






3,3,9 


10,3 


4,7,3 


8,5,8 rat. 




8,5 






6,3,9 


10,3 






4,4,9 


9,4 






3,3,10 


10,3 



Table 1: Examples of singularities and their universality classes of three-matrix 
models with (li,l 2 ,h) = (4, n, 3) 

4 Four-matrix models and beyond 

If the model involves four or more matrices we must always start from an ansatz 
with free parameters such as (four matrices) 

Z n=0 

m 2 = (l 3 -l)(k-l) + (h-l)-X 2 (4.2) 
The two "internal" Schwinger-Dyson equations 

£?! + £?3 = V 2 \B 2 ) 

B 2 + B 4 = Vi(B 3 ) (4.3) 

admit two compensation mechanisms of type I or one of type I and the other of 
type II as illustrated in the four graphs: 
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A 



A 



A4 — 2A3 



1 2 3 4 i 1 2 3 4 i 



A 



A 



A3 — 2A2 — 2A4 



A2 — A3 



1 2 3 4 i 1 2 3 4 i 



Without doubt the universality classes are always [p, q] , this conclusion is 
drawn solely from the possible compensation types. 

In the five-matrix-model case the conclusion is the same, but the arguments 
are more general. Namely, there are up to three types of compensations of type 
II possible (see the graph) 



A 



3A 3 = 6A 2 = 6A4 = 2A 5 



Of course, for a compensation one needs simply Aj = nAj±i, n e IN, the explicit 
cases given here (as A 4 = 2A 3 ) are just examples. 
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One can imagine that what was shown with these examples holds true for 
all /-matrix models with arbitrary /. Thus the universality classes are always 
of the two-matrix model type. The actual calculation proceeds by solving the 
first Schwinger-Dyson equation ( |4.3| ) with the ansatz (3~T) and the method of the 



three-matrix model. Then {pm\ Pm} an d {g^} result as polynomials in {a^-*}. 
Next we have to adjust the {a^} and the {gff} so that 

1. 

r^\z) = V 3 '(r {3 \z)) ~r {2 \z) 

p { £ = for - (l 3 - l)(l A ~ 1) < m < -2 (4.4) 

2. r^'(z) has a maximal zero of degree A4 at z = 1. 

For simple four-matrix models this program can even be performed by hand, e.g. 
for (k,k,k,k) = (3,3,3,3). 
In this case we find a 

[Ai,A 2 ,A 3 ,A 4 ] = [5,4,5,4] (4.5) 

with 

/12 = 2, g>3 — 1 (by normalization) (4-6) 

/i 3 = l, g { ? = ^3(17^-2477-5) 

gf = -L(-^-6»7 + 15) (4.7) 
877^ 

where r] is any of the three real solutions of 

r/ 3 + 7r/ 2 - 15r/ - 5 = (4.8) 

Examples how the different compensations, that we presented in this section ac- 
tually occur, can be seen in appendix B. 



5 Summary 

An algorithm has been described by which all critical points of /-matrix models 
with polynomial potentials in the double scaling limit N — > 00 can be derived. 
These critical points form universality classes [p, q] where p and q are the orders 
of differential operators of a generalized KdV hierarchy and satisfy only 

P > q, p/q ^ IN- 
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We have given a new constructive argument to move this. 

The algorithm treats the Dyson- Schwinger equations, derived from the or- 
thogonal polynomial approach, perturbatively by expansion in power series of a 
rational power of For / > 3 these Dyson-Schwinger equations split into (/ — 2) 
"internal" and (two) "external" equations. A leading order analysis of the internal 
equations fixes the critical point. The external equations determine only the two 
external critical potentials. The fact that only two differential operators result 
from / matrices {Bi}{ =l in the double scaling limit is the result of / — 2 compen- 
sations. We distinguish between two different compensation mechanisms (type 
I and type II) that act on any triplet B i} B i+ %, B i+2 . The differential operators 
of order p and q arise after the compensations and at higer perturbative order. 
Thus we must be able to push the perturbative order to any desired value. This 
fact was responsible already for the derivation of all critical points of the two- 
matrix models in |J (including some not seen before) and is even more decisive 
for the analysis of multi-matrix models. As usual the perturbative order of the 
commutator 

[51,52] 

is p + q — 1 and gives the string susceptibility exponent 7 by 

-2 

7 = • 

p + q — 1 

In the case of the three-matrix models more detailed results have been ob- 
tained. We determined the maximal (i.e. parameter free) critical points and 
found classes of "rational" and "algebraic" critical points, referring to the values 
of their critical coupling constants. The rational class derived from (|3.5| ) is absent 
however, if 

e IN or E IN 

{3 — 1 *i — 1 

where Zi, I3 denote the degree of the external critical potentials. 
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A A linear algorithm 

We make the ansatz 



(z - \) h 



r W (*)= V * ;U (A.1) 



" ' ' \m-h-\, 



We have the "internal" Schwinger-Dyson equation 

r {1 \z) + r^(z) = V^ 2 \z)) (A.3) 

Then ( |A.3| ) amounts to 

f pff, -(i2-l)(i8-l)<m<-2 

E 0i 2) 4n*-l = P« + P^, -1 < m < +1 (A.4) 

*=i I p£>, 2 <m< ^-1)^-1) 

with 

A,^(-i)--)-»g_ + ;- + g (a,) 

In order to produce zeros of order Ai, A3 at z = 1, the coefficients p$,p$ must 
fulfil 



E ]P^=0, 0<n<A!-l 

m=l \ n / 
(Ji-l)(ia-l) 



E ?l°m=0, 0<n<A 3 -l (A.6) 



m=-l \ n 



where we want to maximize Ai, A3. 
With the shorthands 



£' 3) (A.7) 



the m-summations of ( A.2|) can be done analytically to yield 



h-i 



E 9^n k = J2?> 
,(2) 



E C+l*n fe = (A.8) 



fc=0 
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with 



n,k 



■1) 



\ K* , 



+ ( _ 1) A' 1+ l^+^3 



+ (-1) 



V K 3 + 1 , 

Jti+i+r. (K x + K 3 -n 
{ K 3 -l 



(A.9) 



ri.fc 



Kl+n+1 K 1 + K 3 \ 
\Ks-l) 
K x + K 3 -n- T 

^3-1 



("I) 

+ (_!)*!+« 



(A.10) 



and 



We have (a e {1,3}) 



^ = k(h-l) 
K 3 = k(l 3 -l). 



(A.ll) 



(A.12) 



and 



+ $ n _ 1<k ) = (-l) n+1 (^ n>fc + ^ n _! 
= ("!) 



Kl , ATx + K 3 - n - V 



K 3 -2 



(A.13) 



The critical coupling constants {g^} are then determined from homogeneous 
equations 



h-i 



+ $»-i,fc)02i = °» ™>3 



(A.14) 



fc=0 



with ( |A.13|) inserted. All p^\p^ can then be calculated from ( |A.4|) except say 



_(3) _(3) (3) 
P+liHO iP-l 



But these follow from the first equations (|A.8|) for ti G {0,1,2} since R^i\ R 
determine these parameters uniquely. 
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Now we consider the case (l± — 1) = m{l 3 — 1), m G IN\{0, 1}. Then we have 
Ki = mK 3 and the binomial coefficients in ( |A.13| ) have the form 



-ir($ n , fc + $ 



n—l,k ) 



\mK 3 



\m + l)K 3 - n - f 



(A.15) 



and especially $ 2 ,fc can be expressed in terms of (§ n ,k + &n-i,k), n > 3: 



JTi+l 

n=3 



(m + l)if 3 - n - 1 N 
^3-2 



where a 3 = a 4 



1 and a 



m+l 



/(m+l)A- 3 -3\ 

rl if,-i J 

m + l. From this we see that 



(A.16) 



4 3) = _>eW 



(A.17) 



k=0 



can be written as 



m+l Z2-I 

4 3) = E «>(-i) m E + ^-i,*) 



(A.18) 



i=3 



fc=0 



Because of ( |A.14| ) the r.h.s. of (|A.18 ) is zero and R2 is zero, too. This means 

/o\ ______ 

that p_{ is zero (see ( |A.7| )) which is not allowed because then the "external" 
Dyson- Schwinger equations can not be solved. 
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B Maximal critical points for the (4,3,3) model 



To illustrate the fact, that for the f-matrix models with / > 3 one can have more 
than one maximal critical point, we discuss the three maximal critical points 
([4,5,4], [3,3,6], [4,4,5]) of the (4,3,3) three-matrix model. 



The critical point [4,5,4]: 

The functions r^\z), i = 1,2,3 are 

112128152 /T ^ X 
riW = 15?-33 + F-3l + 3-5^ (B1) 

r 2 (z) = -\ + --10 + 10z-5z 2 + z 3 (B.2) 

z l z 

r 3 (z) = --- + 7 --^ z + 9z 2 -7z 3 + 3z A --z 5 + -z e (B.3) 

iK ' 3z 3 5 3 15 K ' 

and we have the critical potentials 

3355 25 2 125 , 125 , ,„ . 

VAx) = x x 2 + x 3 x 4 B.4 

n ; 216 48 8 32 v ; 

V 2 (x) = lx 2 + ^x 3 (B.5) 

72 27 

V 3 (x) = —x + —x 2 -3x 3 . (B.6) 
5 2 

In this case the operators R± and — R3 are equal and therefore the two commu- 
tators [R 2 ,Ri], [i?3,-R 2 ] are obviously equal and this critical point leads to the 
universality class [5,4]. 



The critical point [3,3,6]: 

The functions r^(z), % — 1, 2, 3 are 

la 13 a la 77 a 143 17 

riW = -27 + t?-io?-io7 + ir*-T^ (R7) 

12 a 26 a 11 61 9 , 5 , ,„ , 

roiz) = 1 a H az az H — az (B.8) 

5 z 2 5 z 30 10 2 6 

125 a 18625 " 8125 23125 2 14375 , 

rr>(z) = a H az az H az 

31 ; 432 z 10368 1728 3456 2592 

9125 4 125 5 625 6 

; «z 4 H az 5 a/ (B.9) 



3456 192 10368 
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and we have the critical potentials 

1948250 267493 2 1 6400 6400 3 

1[X) ~ 397953 aX+ 353736 X + 3 ( 14739 + 14739 a)x 

1 , 1250 625 . a .„ , . 

+- a )x 4 B.10 

4 V 14739 4913 ; V ; 
1 19^ 9^ 

V 2 (x) = -x 2 + -(-—a - —)x 3 (B.ll) 

v J 2 3 V 144 144 ; v ; 

Tr/ , 1393 26568 2 1,4478976 4478976, , .„ . 

V 3 (x) = ax-\ x 2 + -( a-\ )x\ (B.12) 

3K J 80 625 3 V 78125 78125 ; ' V ; 

with a being a solution of 

2a 2 + 2a-l = 0. (B.13) 

This is one of the nonrational cases. If one proceeds to higher order one can get 
higher order algebraic equations for more than one variable in order to establish 
a certain critical point. This means that calculations become very complicated or 
even impossible. The operators in this case can be written as 

R x = R 13 + a -7 #1,4 + a" 27 R 1;5 + a" 37 R lfi + a~ 47 R 1)7 (B.14) 
R 2 = R 2t3 + a-m 2A + a- 2 ^R 2 ^ + a- 3 m 2t6 + a-^R 2 j (B.15) 
R3 = #3,6 + 0-^3,7 (B.16) 

where the second index denotes the highest order of the differential operators. 
Because of R\ 3 = R\ 3 = R 3 $ both commutators vanish to leading order. But in 
higher orders one finds 

[R 2 ,R,] = -a- 47 [(-^a-if){# 2 ,3,# 2 ,4} + #i,7-#2,7,#2, 3 ] (B.17) 
[R 3 ,R 2 ] = a-y[{R 2>3 ,R 2A } + R 3! r,R 2>3 }. (B.18) 

During the calculation of these commutators we have used the Dyson-Schwinger 
equations which provide us with equalities such as #1,4 = #2,4,-^1,5 = #2,5, etc.. 
The two commutators ( |B.17| , B.18|) are then commutators of an operator of order 



three with one of order seven and they are equal to the operators of the (7,3) two- 
matrix model, as expected. This critical point then gives a solution that belongs 
to the universality class [7,3]. 



The critical point [4,4,5]: 
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The functions r^\z), i = 1, 2, 3 are 



! + §/?-§/?* 2/3 2 - 4/3 -4/3 2 -2 + 10/3 
n (z) = ^ — 5 + — + 

z 4 Z 6 z l 

+ 4/3 2 + 4 - 12/3 _ 2 _ 3 + + 2 2 + 4 _ 8 

z 5 5 5 

r 2 (z) = -4 + ^^-6/3-4 + (4/3 + 6>+(-4-/3)^ 2 + z 3 (B.20) 

-11+4/3 2 „ rt „ rt „„, , 926 28 ^ 362 
r 3 (z) = ^— + 69 + /3 - 26/3 2 + ( — -— (3 + — (3 2 )z 

Z OOO 

+(2 76 + 13(3 - 112(3 2 )z 2 + (-247 - 16(3 + 104(3 2 )z 3 
+(11(3 - 58(3 2 + 133)/ + (-4(3 + 18(3 2 - 40)/ 

+(^-f^ 2 + f)/ (B.21) 

and we have the critical potentials 

Tr , x 11937^ 294725 263361 n 1 , 8013 „ 2 157853 1221 „. 2 

Vi(x) = /TaH x-\ (3x--( (3 2 1 (3)x 2 

w 19652^ 39304 39304 1 2 V 19652 M 39304 19652 M; 

1,22959 „, 154635 161649 ^ , 1,7933 27915 „ 4254 4 ,„ 

( R 1 3)x A ( 3 3 )x (B 22) 

3 V 39304 M 19652 39304 f J 4 V 9826 39304^ 4913 M ; V ' ; 

V te) = \* + \(\P-T? + ^ (R23) 

Tr/ . 115848 „ 378696 2134863^ 1 ,21225 219258^ 62334 2 

V 3 (x) = 3 2 x-\ x-\ (3x+-( 1 3-\ 3 )x 2 

31 ; 222605^ 222605 222605 ^ 2 1 44521 44521 ^ 44521 ^ ; 

1, 34521 n 17040^ 4568 . , ,„ 

+ 3<-4452l' 3 -4452l' 3 -44521^ (&24) 
with (3 being a solution of 

3(3 3 -7(3-1 = 0. (B.25) 

The operators in this case are 

R-y = R 1A + a- J R 1>5 (B.26) 
R 2 = R 2:4 + a-' r R 2j5 (B.27) 
Rs = ^3,5 (B.28) 

The leading order of the first commutator vanishes again (i?i,4 = #2,4) but the 
second commmutator does not: 

[i^fli] = -a- J [R 1 , 5 -R 2 , 5 ,R2,4] (B.29) 
[^3,^2] = [#3,5,^2,4]- (B.30) 

Then we have — i?i ;5 + i? 2 ,5 = -^3,5 as operators of order five and find the univer- 
sality class [5,4]. 
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C Examples for the three-matrix model 



In this appendix we give the potentials that belong to rational critical points of 
the type that is discussed in Appendix A. This means that A2 assumes its maximal 
possible value {l\ + 1 3 — 2). Models of the class (5, x, 3) were omitted because they 



belong to the cases where the construction based on the ansatz (|3.5| ) fails. 



model 


critical point 


critical potentials 


3 2 3 


[3 4 31 


Vi = —3 T - T 2 + - T 3 

V <_J iXj JU I g >x> 

Vz = V l 


3, 3, 3 


[4 4 41 


T/ HI t _i_ 27 ™2 _i_ 3 ™3 

1 16 8 ~ r 4 ,x ' 

^3 = ^1 


3 4 3 


[5 4 51 


T7" _ 1305 , 75 2 , 25 3 
vi 112 ~ 16 ' 48 
y _ 1 ~3 3 4 
v 2 3 140 X 

^3 = ^1 


3, 5, 3 


[7 4 71 


_ 51485 6125 2 , 1225 3 
1 3072 512 x "r 768 x 
\r _ 1 ™3 1 ™4 _i_ 1 ™5 
3 10 X ^ 210 X 

14 = ^i 


3, 6, 3 


[8 4 81 


t /" 957285 , 9261 2 i 147 3 

Vl 45056 ~ r 2048 * r 1024 X 

I/„ _ 1 T 3 3 4,9 5 1 6 
v 2 — 3 x 4 x -r 7Q x 19g x 

^3 = ^1 


4 2 4 


f3 6 31 


vi 27 18 x ~ 27 x ~ 108 x 
^2 = |X 2 

^3 = ^1 


4, 3,4 


[5, 6, 5] 


T/ _ 6595 125 2 , 250 3 125 4 
vi lag ^ 126 x ~ 27 108 

14 = ^1 


4, 4,4 


[6, 6, 6] 


t /" _ 311127320 , 287800 2 i 160000 3 i 2000 4 

1 5845851 1 59049 1 59049 1 19683 
^2 = |x 2 + A a; 3 + _I_ a; 4 

^3 = ^1 


4, 5,4 


[7, 6, 7] 


t r 508518850 11502260 2 , 266200 3 2662 4 

Vl 6908733 X 767637 X 1 59049 X 19683 X 

\f 1 rp'i 1 5 ™4 1 1 „5 

v 2 3 264 4862 ^ 
^3 = ^i 
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4, 6,4 


[9, 6, 9] 


y — 300701859304 , 3004141294 2 , 204526784 3 i 7304528 4 
Vl 3171108447 1 62178597 1 14348907 1 14348907 
T/ — 1 ~3 1 1 ~4 1 3 „5 1 1 ~6 
V2 — 3 X 16 ^ 1430 X ^ 54264 X 

V 3 = Vi 

* o * J- 


5, 2,5 


[3, 8, 3] 


t/ _ 24567 103 2 3 3 3 4,1 5 
1 400 200 1000 2000 1 50000 

\/ 2 = |o; 2 

V 3 = Vi 


5, 3, 5 


f5, 8, 51 

L ? ? J 


t /" _ 109122293 5831 2 , 232211 3 , 21609 4 , 2401 5 
1 720896 16384 1 24576 1 8192 1 20480 

V2 2 x 4 6 2 x 

V 3 = Vi 


5, 4,5 


[7, 8, 7] 


t /" 5926366827535 500520251 2 i 50463788221 3 
Vl 25518145536 327155712 1 981467136 
175765205 4 i 35153041 5 
16777216 1 83886080 
T/ _ 1 T 2 1 3 , 1 4 
y 2 2 63 ~ 19448 

v 3 = vi 


5, 5, 5 


[8, 8, 8] 


t /" 318301318564336604805 , 195795357608389 2 
Kl 1020323700831944704 1 30374008717312 X 
, 45439661273853 3 i 3691069305 4 , 35153041 5 


1 7593502179328 1 17179869184 1 21474836480 
U _ 1 T 2 3 3,3 4 1 5 
V2 ~~ 2 X 14 X ^ 1144 X 148580 X 

V 3 = Vi 


5, 6, 5 


f9, 8, 91 

Y 1 ' J 


Y - 489711657601739357169 , 2079610272730221 2 
1 1235128690480775168 1 60748017434624 
, 1697130639464183 3 i 93493458289 4 i 1908029761 5 


1 91122026151936 1 137438953472 1 343597383680 
T/ _ 1 ™3 7 4 1 7 5 7 6 
^ 2 3 572 ^ ~ 74290 ^ 35357670 x 

v 3 = V 


4, 2,3 


[3, 5, 3] 


Vi = 8 x - x 2 + | x 3 + \ x A 
V 2 = \x 2 

TA — T - T 2 — /t-3 

3 4 2 12 


4, 3,3 


[4, 5, 4] 


t/ _ 3355 _ 25 _2 , 125 „3 125 _4 
Vl — 21g X 48 8 32 X 

V 2 = \x 2 + ^ 

t / 72 , 27 2 q 3 

1/3 = yi + yr-or 


4, 4,3 


[5, 5, 5] 


T /" 223600 , 1940 2 , 8000 3 1 250 4 

f 1 926I x "r 343 1029 343 

V2 2 15 120 *^ 

V 3 = f x-25x 2 -f|a; 3 


4, 5,3 


[6, 5, 6] 


t /" 4014340 154396 2 , 53240 3 2662 4 

1 120393 4459 1 1029 343 

V 2 = 3 x 3 + 24 a; 4 + 100T ^ 5 
V 3 = -56x+^fx 2 -f x 3 
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4, 6,3 


[7, 5, 7] 


t /" 43138744 , 1969880 2 , 2725888 3 , 1362944 4 

Kl 1008423 1 19773 1 19773 1 59319 
T/ — 1 ~3 , 9 ~4 , 3 ~,5 , 1 ™6 
v 2 3 X "r 64 x " r 286 X ^ 4896 X 

\/ 3 = _9|Za; _ 324x 2 -48x 3 


5, 2,4 


[3, 7, 3] 


T/ _ 8331 21 2 , 3 3 3 4, 1 5 
1 256 32 ~ 128 128 ~ 1280 

T/ o _ 626 11 2 1 3 , 1 T 4 
v 3 27 24 27 ~ 864 


5, 3, 4 


f4, 7, 41 


t/ _ 694729 1029 2 , 28126 3 , 21609 4 , 2401 5 
1 10000 5000 1 1875 1 2500 ~ l ~ 3125 

V2 2 *^ 210 

t / 26045 1375 2 125 3 125 4 

V 3 X g X 


5, 4,4 


[5, 7, 5] 


t/ 239105243531 4369673 2 1 10962726929 3 
Vl ' 2190240000 780000 1 63180000 
35153041 4 1 35153041 5 
405000 1 4050000 
T/ _ 1 ~2 2 3, 1 „4 
V 2 — 2 X 63 X ^ 4576 X 
\/ 3 = 12680 + 2144 ^2 512 ^3 + 16 ^4 
J 77 7 3 


5, 5,4 


[6, 7, 6] 


t/ 140287221613 , 665951209 2 , 80349808 3 
1 940970706 1 8739666 1 336141 X 

, 175765205 4 , 35153041 5 

1 3084588 1 11567205 
y _ 1 T 2 3 3,3 4 1 5 
^ 2 2 14 572 35530 
t / 326419 14641 2 6655 3 1331 4 
V3 37g X 14 X 27 x 108 X 


5, 6, 4 


f7, 7, 71 

L ' > J 


y _ 2615115300936941 7713352643381 2 , 98261624661739 3 
v l 14045490038528 85875958116 1 2748030659712 

4581179456161 4 , 4581179456161 5 

2748030659712 X 1 247322759374080 
T/ _ 1 T 2 16 3 , 45 4 4 5 , 1 T 6 
y 2 2 7 ~ r 286 1615 ~ 91770 
t/ 1220198 | 62209 2 343 3 , 343 4 
V 3 — 297 X 1 72 X 9 X " r 864 ^ 


6, 2, 3 


f3, 7, 31 


Vi = 24x - x 2 + 3 x 4 + 1 x 5 + 1 x 6 

1 '5 '6 

V2 = |X 2 


6, 3, 3 


[4, 7, 4] 


t /" _ 4217031 729 2 , 40743 3 308367 4 
1 100000 1600 1 400 320 
, 177147 5 19683 6 
1 160 64 

V2 = 2 X 2 + 135 

\/ 3 = -^f x + f x 2 -3x 3 


6, 4,3 


[5, 7, 5] 


T/ _ 1157614349283 6163075107 2 , 20496602637 3 
1 18564650000 11881376000 1 14851720 
1 4321777591467 4 , 80387359983 5 , 80387359983 6 
1 475255040 1 7425860 1 23762752 X 

Vo - i X 2 + — X 3 + — X 4 
r 2 2 X ^ 81 X ^ 1680 X 

V 3 - 44 2 5 8 ° 5;r 675 ^ 
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6, 5, 3 


[6, 7, 6] 


Y — 34173103071667 18059440511 2 i 3311338930137 3 
Vl 408422300000 32673784000 1 1633689200 
684760877857 4 i 41615795893 5 41615795893 6 
118813760 1 11881376 71288256 

\/ 3 = -10695 x + ^x 2 -^^ 


6, 6, 3 


[7, 7, 7] 


t /" 788322756540293184 , 3521389972070481 2 , 865581376911776 3 

1 7468028442184375 1 298721137687375 1 2757425886345 
, 1969585683008454 4 , 242703321647976 5 1 3370879467333 6 
1 11948845507495 1 11948845507495 1 4779538202998 

T/ — 1 „2 1 32 3 1 3 4 1 8 5 1 1 ™6 

r 2 2 x-r g x-r 7 x-r 665 x-r 10530 x 

t /" _ 180675 11907 2 1323 3 
v 3 4 x 4 x 16 


6, 2,4 


[3, 8, 3] 


Y — 190536 851 2 16 3 1 1 -4 
V1 3125 1250 1875 ^ 50 

1 8 ~,5 1 1 ™6 

1 3125 1 18750 

T/ — 174 „ 3 2 , 1 „3 , 1 „4 
5 x s X ^ 75 X ^ 4000 X 


6, 3,4 


[4, 8, 4] 


t/ 775352871324 3256929 2 , 41065704 3 52337097 4 
Vl 6355146875 20336470 X 1 4621925 X 1848770 X 

1 708588 5 19683 6 

1 84035 33614 

V2 = 2 x2 + 1485 ^ 

t/ 5825 „ 4325 „2 , 125 „3 125 „4 
v 3 — 27 24 ~ l ~ 2 32 


6, 4,4 


[5, 8, 5] 


Y 45952369797 986577 2 1 2048274 3 1 859977 4 
Kl 244268750 1645600 1 4675 1 880 
+ 1944 5 + 81 6 

T/ _ 1 T 2 , 8 3 , 1 4 
f 2 2 ^ ~ 495 ~ 15504 

y 3 - 169600 x + 7550 3.2 + 8000 ^3 + 25Q % A 


6, 5,4 


[6, 8, 6] 


Y — 2294017899932829 30776961716517 2 , 1636161080121 3 
Kl 8980595525000 1436895284000 X 1 780921350 X 

320385186351 4 , 1203384114 5 200564019 6 
113588560 1419857 2839714 
T/ - _ 1 r 2 1 4 3 1 1 „4 1 6 5 
^2 2 ^45 ^816 ^ 1562275 

Y — 8025008200 227667550 2 , 13310000 3 332750 4 
3 120393 4459 1 1029 343 


6, 6,4 


[7, 8, 7] 


t /" _ 154075511895816746448 , 23090824363335471 2 , 9611096332631616 3 
1 476632759170528125 1 131484899081525 1 5716734742675 
1 928245980327382 4 , 112290178445568 5 , 779792905872 6 
1 1143346948535 1 1143346948535 1 228669389707 

T/ _ 1 ™2 , 16 3 I 3 4 , 24 5 I 1 ™6 

V2 2 X 27 x 136 x 120175 x 1 1941840 

V 3 - 100 f 7 806 x + 254 2 6 f 3 x 2 + 10648 x 3 + 13 f x 4 
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